% This is a script for reading in the phantom files, as an example
Nx = 495; Ny = 615; Nz = 752;

% Load in optical and acoustic phantoms
fid = fopen('./MergedPhantom.DAT', 'r');
phan = fread(fid, 'uint8=>uint8'); phan = reshape(phan, [Nx, Ny, Nz]);
fclose(fid);



